Cycloaddition Functionalizations to Preserve or Control 
the Conductance of Carbon Nanotubes 
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We identify a class of covalent functionalizations that preserves or controls the conductance of 
single-walled metallic carbon nanotubes. [2+1] cycloadditions can induce bond cleaving between 
adjacent sidewall carbons, recovering in the process the sp 2 hybridization and the ideal conduc- 
tance of the pristine tubes. This is radically at variance with the damage permanently induced by 
other common ligands, where a single covalent bond is formed with a sidewall carbon. Chirality, 
curvature, and chemistry determine bond cleaving, and in turn the electrical transport properties 
of a functionalized tube. A well-defined range of diameters can be found for which certain addends 
exhibit a bistable state, where the opening or closing of the sidewall bond, accompanied by a switch 
in the conductance, could be directed with chemical, optical or thermal means. 



Chemical functionalizations of carbon nanotubes 
(CNTs) are the subject of intensive research PIQ, and 
could offer new and promising avenues to process and as- 
semble tubes, add sensing capabilities, or tune their elec- 
tronic properties (e.g., doping levels, Schottky barriers, 
work functions, and electron- phonon couplings). How- 
ever, the benefits of functionalizations are compromised 
by the damage to the conduction channels that follows 
sp 3 rehybridization of the sidewall carbons |3,|4|, as ev- 
idenced by absorption spectra and electrical transport 
measurements 0, 0, ■ We report here on a class of cy- 
cloaddition functionalizations that preserves instead the 
remarkable transport properties of metallic CNTs. In ad- 
dition, we identify a subclass of addends that displays a 
reversible valence tautomerism that can directly control 
the conductance. 

We focus here on [2+1] cycloaddition reactions, where 
the addition of a carbene or a nitrene group saturates 
a double-bond between two carbon atoms, forming a 
cyclopropane-like three-membered ring. Such function- 
alizations have been reported extensively in the litera- 
ture 0, 0, ^3 • All our calculations are performed using 
density-functional theory in the Perdew-Burke-Ernzerhof 
generalized gradient approximation (PBE-GGA) ^lj, ul- 
trasoft pseudopotentials [T^j, a planewave basis set with 
a cutoff of 30 Ry for the wavefunctions and 240 Ry for 
the charge density, as implemented in Quantum-ESPRESSO 

We examine first the simplest members in this class of 
addends, CH2 and NH. Fig. da) and (b) show the two 
inequivalent choices available on a (5,5) metallic CNT; 
for convenience, we label these as "S" (skewed) and "0" 
(orthogonal) , reminiscent of the relative positions of the 
sidewall carbons with respect to the tube axis. Our sim- 
ulation cells include 12 n carbons for a given (n,n) CNT, 
plus one functional group. We use a 1x1x4 Monkhorst- 
Pack mesh (including T) for structural optimizations, and 
a 1x1x8 mesh for single-point energy calculations, with 
a cold smearing of 0.03 Ry [14j. 



First, and for (n,n) metallic tubes, we highlight how 
strongly the reaction energies of these functionalizations 
depend on the curvature of the nanotubes, and on their 
attachment sites S and 0. We plot in Fig. |2ta) the reac- 
tion energies AScnt (defined as AEcnt — -EcNT-func — 
Egnt — -Efunc, func=CH2 or NH), taking as a zero refer- 
ence the same reaction on graphene. The reaction ener- 
gies have a well-defined linear dependence on curvature, 
clearly demonstrating the higher reactivity of smaller- 
diameter tubes. The site is always more stable, and 
significantly so for all diameters considered; at room tem- 
perature all small-diameter armchair CNTs strongly fa- 
vor the configuration. The energy difference between 
the and the S form of the (5,5) CNT is 1.24 eV, which 
is in good agreement with other planewave basis calcula- 
tions (1.24 eVjRef. 0]) and localized basis calculations 
(1.4 eV, Ref. 0). 

Second, we find that the Ci-Cg distance for this con- 
figuration (die m Fig- E^ c )) is much larger than the 
usual C-C distance (1.54 A in diamond and 1.42 A in 
graphite), a clear indication that the sidewall bond is 
broken [H [H El . For the CH 2 and NH cycloaddi- 
tions, we observe bond cleaving for all nanotubes studied 
(up to (12,12)); on the other hand, in graphene the bond 
is intact. We estimate the critical diameter that sepa- 
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FIG. 1: The three different configurations for a func- 
tional group on an armchair nanotube are shown (CH2 on 
a (5,5) CNT): (a) skewed S, (b) orthogonal with an intact 
("closed") sidewall bond, and (c) orthogonal with a broken 
("open") sidewall bond. 
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FIG. 2: (a) Energy change AEcnt upon functionalization 
as a function of curvature, for (n,n) armchair CNTs - the zero 
reference is for graphene. (b) Sidewall equilibrium bond dis- 
tance C1-C6 (die) for (n,n) CNTs and for a bent graphene 
sheet, functionalized at the stable site. The sidewall bond 
is broken in all the (n,n) CNTs considered. Continuous bend- 
ing of a graphene sheet shows that a well-defined transition 
from the closed to the open form takes place as the curvature 
increases. 



rates the two regimes by bending a graphene sheet: we 
can see in Fig. |^b) that such model closely reproduces 
the nanotube results, and a sharp transition from the 
bond-intact to the bond-broken form takes place around 
a diameter of 2.4 nm, i.e. a (18,18) tube, for CH2 func- 
tionalizations. 

Broken or intact sidewall bonds play a fundamental role 
in the electronic transport properties of a metallic nan- 
otube. We show this in Fig. |3 where the Landauer 
conductance is calculated for a (5,5) CNT functionalized 
with dichlorocarbene (CCI2) (first reported experimen- 
tally in 1998 0), using an approach recently introduced 
by us that allows to treat realistic nanostructures with 
thousands of atoms while preserving full first-principles 
accuracy For this diameter the sidewall bond is bro- 
ken, and the scattering by a single CCI2 group is found to 
be remarkably weak, with the conductance approaching 
its ideal value. Even after adding 30 groups on the cen- 
tral 43 nm segment of an infinite nanotube the conduc- 
tance is only reduced by 25%. This is in sharp contrast 
with the case of a hydrogenated tube, where the conduc- 
tance drastically drops practically to zero when function- 
alized with a comparable number of ligands. This result 
is easily rationalized. Hydrogen and other single-bond 
covalent ligands induce sp 3 hybridization of the sidewall 
carbons, and these chemical defects act as very strong 
scatterers Q . Such dramatic decrease in the conductance 
has also been recently confirmed experimentally 0. On 
the other hand, after bond cleavage Ci and recover 
a graphite-like bonding environment (Fig. ^c)), with 
three covalent bonds to their nearest neighbors. Their 
electronic orbitals go back to sp 2 hybridization, allowing 
for the p z orbitals of Ci and Cq to be recovered (as con- 




FIG. 3: Quantum conductance for a (5,5) CNT functional- 
ized with one CCI2 group, with 30 CCI2 groups, and with 30 
hydrogen pairs (dashed line: pristine (5,5) CNT). A CCb ad- 
dend will choose the O-open configuration, while hydrogens 
prefer to pair in the S configuration; the results corresponds 
to these stable choices. The conductance for the energetically- 
unstable O-closed configuration is also plotted. The 30 func- 
tional groups are arranged randomly on the central 43-nm 
segment of an otherwise infinite tube; the conductance is then 
averaged over 10 different configurations. 



firmed by inspection of the maximally-localized Wannier 
functions |l9l |) and to contribute again to the graphitic ir 
manifold. The net result is that conductance approaches 
again that of a pristine CNT, highlighting the promise of 
cycloadditions in preserving the conductance of metallic 
nanotubes. 

In principle, around the critical diameter shown in Fig. 
|2Ib), a functional group could be found that displays 
as stable states on the same tube both the open and the 
closed configuration. If the case, interconversion between 
the two valence tautomers would have a direct effect on 
the conductance. We illustrate this paradigm with the 
dichlorocarbene example of Fig. [31 if we force the tube 
in its closed configuration, with all the sidewall bonds 
frozen in their ideal pristine-tube geometry, conductance 
decreases by a factor of 2 or 3 with respect to the case 
of the relaxed tube [20j. The configuration where the 
sidewall bonds are intact is not stable for this case, but, 
depending on the chemistry of the addends and the di- 
ameter considered, optimal ligands could be found for 
which a double- well stability is present. 
In order to identify the factors that determine stabil- 
ity in the open or closed form, we first screen several 
addends on small molecules. We find that the bridged 
l,6-X-[10]annulene (inset of Fig. 0fa)) is an excellent 
molecular homologue of a functionalized CNT |17|. It is 
well-known that the substitutional group X dictates the 
preference for the annulene (henceforth labelled as lo) 
or for its valence tautomer, a bisnorcaradiene derivative 
(lc), corresponding to the open and closed configura- 
tions of a functionalized CNT [H |H |H |H Hi H3- 
Similar tautomcrization between an open (2o) and a 
closed (2c) form takes place in a pyrene derivative 
[H (inset of Fig. gfb)). We thus tested on these 
molecules the substitutional groups X= CH2, NH, SiH2, 
C(N0 2 ) 2 , C(CN) 2 , C(CCH) 2 , C(CH 3 ) 2 , C(COOH) 2 , 
CC1 2 , C(NH 2 ) 2 , C 6 H 4 29], and d 3 H 8 Q. To assess 
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FIG. 4: Potential energy surface as a function of die for select 
cases of (a) 1 and (b) 2. C(CN) 2 (violet) and C(CCH) 2 (blue) 
show a double-well minimum in both 1 and 2. C(N02)2-rot 
(dashed red) indicates the unstable conformation where the 
two NO2 groups are rotated from their equilibrium position 
(solid red) by 90° . 

the accuracy of our PBE-GGA approach, we compare 
in Table [I] our results for die m 1 with those obtained 
from experiments or other theoretical methods (second 
order M0ller-Plesset perturbation theory (MP2) |30| and 
Becke thr ee p arameter Lee- Yang-Parr hybrid functional 
(B3LYP) |31|), finding excellent agreement. Hydrogen 
and halogens have been reported to stabilize lo both ex- 
perimentally and theoretically 0I1EI13- Cyano group 
favors lc experimentally (the possibility of coexistence 
with lo is discussed) |27|, and two minima have been 
actually predicted theoretically [13, HH . 
We show the potential energy profile of select groups in 
Fig. 0] All carbenes stabilizing the closed form in both 
1 and 2 share a common feature: partially-occupied p 
orbitals parallel to the in-plane p a orbital of the bridge- 
head On atom 36] (the plane considered is that of Ci- 
Cix-C6 of Fig. nfc)). This conclusion is strongly sup- 
ported by examining the energy minimum conformation 
for X=C(NC>2)2- At equilibrium, the two oxygen atoms 
in the NO2 group lie on a line parallel to the Ci-Cg bond, 
and the open form is stable (solid red in Fig. 0J. The 



TABLE I: Experimental and theoretical d\§ of 1. Note 
that the calculations assume isolated molecules at K, while 
experimental data are obtained from crystalline systems at 
finite temperature. Theory predicts two stable minima for 
X=C(CN) 2 . The long die of X=C(CH 3 ) 2 indicates that the 
potential energy surface would be very flat, which is also pre- 
dicted in Fig. 2] especially for 2. 



X 


Expt. 


MP2 B3LYP 


This work 


CH 2 


2.235° 


2.25L 2.279° 


2.278 


CF 2 


2.269° 
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2.300 


C(CN) 2 
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1.572, 2.245 
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NH 
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2.237 s 


2.239 
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FIG. 5: Potential energy surface for (n,n) CNTs functional- 
ized with C(CN) 2 . Both (10,10) and (12,12) CNTs display a 
double- well minimum. 



sidewall bond will switch from open to closed upon rota- 
tion of the two NO2 groups by 90° thereby placing the 
p orbitals of the nitrogens parallel to the p a orbital of 
Cn (dashed red in Fig. 0J. Among all the substituents 
screened, we find that SiH 2 , C(CCH) 2 and C(CN) 2 show 
most clearly the presence of two minima in their potential 
energy surface; we choose here X=C(CN)2 as the most 
promising candidate since both lc 27] and a Ceo deriva- 
tive have already been synthesized. 
We explored, therefore, the potential energy landscape 
for armchair CNTs in the case of C(CN)2 cycloadditions. 
The results shown in Fig. reflect closely those found 
in the molecular homologues. A unique minimum in the 
open form is found in small-diameter tubes, as is gener- 
ally the case in these [2+1] cycloadditions. As the di- 
ameter is increased, the signature of the closed minimum 
starts to appear, first as an inflection point ((5,5) CNT), 
then as a local minimum for the (10,10) CNT (0=1.36 
nm, as in lc), and finally as a global minimum for the 
(12,12) CNT (0=1.63 nm, as in 2c). As discussed be- 
fore, the conductance is controlled by the bonding and 
hybridization of the sidewall carbons. We compare in 
Fig. the two stable open and closed forms for the 
(10,10) CNT functionalized with C(CN) 2 . The scattering 
induced by a single group is negligible, especially in the 
open form, and the conductance around the Fermi energy 
is extremely close to its ideal value (Fig. E{a)). As the 
number of functional groups is increased, the difference 
between the two minima becomes rapidly apparent (Fig. 
lib)). 

Two conclusions can be drawn: First, even with a large 
number of functional groups, the conductance of the tube 
is well preserved, whenever cycloaddition breaks the side- 
wall bond. Second, a subclass of substituents can be 
found (e.g., C(CN)2) that stabilize two tautomeric forms 
on the same tube, separately displaying high and low con- 
ductance. Several mechanisms, including photochemical, 
electrochemical, and thermal, could then direct intercon- 
version between the two tautomeric forms. Photochem- 
ical and electrochemical interconversion rely on the fact 
that the energy levels of the frontier orbitals are affected 
by die, depending on their symmetries and charge distri- 
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dation (DMR-0304019); computational facilities have 
been provided through NSF (DMR-0414849) and PNNL 
(EMSL-UP-9597). 
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FIG. 6: Quantum conductance for a (10,10) CNT functional- 
ized with C(CN)2, in the O-open and O-closed stable configu- 
rations (dashed line: pristine (10,10) CNT). (a) Single group, 
(b) 30 functional groups arranged randomly on the central 
32-nm segment of an infinite tube; the conductance is then 
averaged over 10 different configurations. 



butions (e.g., the bond weakens as a filled orbital that has 
bonding character along Ci-C 6 is emptied) [^,[2^. Both 
photochemical excitations or electrochemical reduction 
or oxidation can populate or depopulate those frontier 
orbitals that favor the open or closed form; as a result, 
they would modulate the bond distance, and ultimately 
the conductance. As a proof of principle, time-dependent 
density functional calculations in 2 for X=C(CN)2 show 
that the first singlet excitation ( Si) drives the system 
from the open to the closed form |38j ■ A similar conclu- 
sion is drawn from experimental observations in 2o for 
X=CH2, where this open form is stable in the ground 
state while the closed 2c is presumed to be stable in 
the Si energy surface [2^. Temperature also plays an 
important role, and 13 C NMR spectroscopy and X-ray 
data have captured the temperature-dependent equilib- 
rium of similar fluxional systems: a higher temperature 
stabilizes lo for X=C(CN)(CH;ih while destabilizing it 

for X=C(CH 3 ) 2 lum. 

In conclusion, our calculations predict that 1) a broad 
class of cycloaddition functionalizations on narrow- 
diameter nanotubes recovers, as a consequence of bond 
cleaving, the conductance of the original pristine tubes, 
allowing for organic chemistry approaches to manipula- 
tion and assembly that preserve the remarkable electronic 
properties of these materials, and 2) that a subclass of 
addends, exemplified in this work by dicyanocarbene, ex- 
hibits fluxional behavior that could be controlled with 
optical or electrochemical means. Such conductance con- 
trol, if realized, could have practical applications in nano- 
and opto-electronics, chemical sensing, and imaging. 
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